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ABSTRACT 

Context. The thick disk rotation-metallicity correlation, <9V^/<9[Fe/H]= 40 4- 50 km s~'dex~' represents an important signature of the 
formation processes of the galactic disk. 

Aims. We use nondissipative numerical simulations to follow the evolution of a Milky Way (MW)-like disk to verify if secular 
dynamical processes can account for this correlation in the old thick disk stellar population. 

Methods. We followed the evolution of an ancient disk population represented by 10 million particles whose chemical abundances 
were assigned by assuming a cosmologically plausible radial metallicity gradient with lower metallicity in the inner regions, as 
expected for the 10-Gyr-old MW. The two cases of a disk with and without a bar were simulated to compare the evolution of their 
kinematics and radial chemical properties. 

Results. Migration processes act in both cases and appear to be enhanced in the presence of a central bar. Essentially, inner disk stars 
move towards the outer regions and populate layers located at higher |z|. In the case of an evolved barred disk, a rotation-metallicity 
correlation appears, which well resembles the behaviour observed in our Galaxy at a galactocentric distance between 8 kpc and 10 
kpc. In particular, we measure a correlation of SV^/dlTe/H]^ 60 km s~'dex~' for particles at 1.5 kpc < |z| < 2.0 kpc that persists up 
to 6 Gyr. 

Conclusions. Our pure W-body models can account for the vs. [Fe/H] correlation observed in the thick disk of our Galaxy, 
suggesting that processes internal to the disk such as heating and radial migration play a role in the formation of this old stellar 
component. In this scenario, the positive rotation-metallicity correlation of the old thick disk population would represent the relic 
signature of an ancient inverse chemical (radial) gradient in the inner Galaxy, which resulted from accretion of primordial gas. 
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1. Introduction 

Two distinct disk populations are observed in the solar neigh- 
bourhood in terms of kinematics, average metallicities, and 
stellar ages. They belong to a thin disk, y ounger than 7-8 Gyr, 
and to a thick disk, older than 8-9 Gyr (e.g. Haywood 200§l and 
references therein). 

We now have several different scenarios for the origin of 
the thick disk: kinematical heating of a pr e -existing old 
disk via minor mergers ( Ouinn et al. I 1 993b Robin et al. I 



1996; IVillalobos&Helmil [2008; O uetalJ 120111) . accretion 
of satellites t hat have deposited their stellar debris in planar 
configuration (lAbadi et al. II2003I). gas accretion a t high redshift 
and stars formed i n situ (iBrook et al. I 120051). wet merger 



a "fossil record" of its evolutionary history. Therefore crucial 
information on the mechanism that dominated the formation of 
the thick disk is encoded in the chemical properties of its stars. 
By comparing the chemical properties of bulge, thick and thin 
disk stars, and by correlating them with other kinematical data, 
one can identify the processes that play a dominant role in the 
formation of the different Galactic components. 
On this matter, an intriguing correlatio n between th i ck disk 
rotation and metallicity was found by ISpagna et al. | ( 20loh . 
Subse quently, this correlation was confirm ed bv | Lee et al. 



(2011 



(2011 



(Matt el et al. I 1201 II) . clumpy disks ( Bournaud & Elmegreen 



20091). and radial m igration (ISchonrich & Binnev I l2009a lbT 



Loebmap et al. 1 1201 ll). T he interested reader can consult the 
work bv lLee et al. I d201 ll) for a detailed overview of the possible 
origins of the thick disk. 

Chemical properties of stars provide important clues to dis- 
entangle the puzzle of the Galaxy formation. By measuring 
abundance ratios of stars in different parts of the Galaxy, one can 
infer how fast metal enrichment proceeded and the timescale 
over which the different regions were formed. The distribution 
of the chemical elements in our Galaxy can be considered as 



) using similar data and, very recently, by Kordon atis et al. 
) using fully independent spectroscopic data. The Spagna 
et al. correlation is rad ically different from what is observed 
for the thin disk stars (lHavwood 1120081) . and from predictions 
made b y the c hemodynamical model s of ISchonrich & Binnev I 
(2009a)' and of lLoebman et al. I (120 1 11) for young stars. 
In the context of stellar migration models, the negative rotation- 
metallicity correlation shown by the thin disk stars in the 
solar neighbourhood derives from the "immigration" of slower 
rotating stars that come from the inner disk regions and are 
characterised by stellar populations with higher chemical 
abundances. 

Conversely, the positive correlation shown by the vs. [Fe/H] 
distribution of thick disk stars observed in situ above 1 kpc 
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from the plane could be explained with a population of much 
older stars that come from the inner regions of the Galaxy and 
have formed during the early disk star formation phases from 
molecular clouds with relatively low chemical abundance and 
higher a-enrichment. Note that relationships between thick 
disks and bulg es in externa l galaxies have been known for quite 
some time (Freemanl ll987l) . 

Although this substantially positive V$ vs. [Fe/H] correlation 
has not yet received much theoretical attention, a weak pos- 
itive correlation is fo und in the simulations carried out by 
lLoebm an et al H d2011l). while t he statistical model developed 
bylSchonrich & Binnev (2009a) indicates only a mild trend (10 
km s _1 dex 1 ) at z =0 kpc, which decreases with height and 
disappears for \z\ > 1 kpc (Schonrich 2009, private communica- 
tion). 

In this paper we show that if we assume the early radial 
chemical gradient derived from the inside-out formation and 
chemical ev olution model of the Galacti c di sk as originally 
sugges ted bv lMatteucci & Francois I dl989h and lChiappini et al.l 
(2001), a positive correlation between rotation velocity and 
metallicity can be established as a result of radial migration 
and heating processes of stars from the inner region of the disk. 
These processes are investigated by evolving a stellar disk with 
N- body simulations and thus following how the chemical prop- 
erties assigned to the initial configuration are redistributed. 
Our model cannot be viewed as a complete galaxy evolution 
model, because we neglected the gradual formation and growth 
of the stellar disk. However, our approach allows us to analyse 
the dynamical evolution of the non dissipative component, and, 
therefore, to separate the chemical evolution from dynamics, i. 
e. leaving any redistribution of chemical gradients to stellar mo- 
tion. 

2. Observed rotation - metallicity correlation 

The observed rotation-metallicity correlation is based on a new 
kinematic catalogue derived by assembling the astrometric pa- 
rameters extracted from the database used for the construc- 
tion o f the Second Guide Star Catalog (GSC-II; lLasker et al. I 
120081) with spectro-photometric data from the Seventh Data 
Release of the Sloan Digital Sky Survey (S PSS DR7; see e.g. 
lAbazaiian et all 120091: lYannvetal.l [2009). The SDSS-GSC- 
II catalogue contains positions, proper motions, classification, 
and ugriz photometry for 77 million sources d own to r ~ 20, 
over 9000 square-degrees. Spasn a etal.1 (|2010) analysed 27 000 
FGK (sub)dwarfs in the solar neighbourhood and found evi- 
dence of a rotation-metallicity correlation, cJV^/SfFe/H]* 40 -r 
50 km s dex -1 , amongst in situ thick disk stars, located at 1 kpc 
< \z\ < 3 kpc and selected with abundance —1.0 <[Fe/H]< -0.5. 
The thin disk contamination is expected to be negligible in this 
sample; however, this is not true for the halo stars, consequently, 
the thick disk and halo populations had to be deconvolved by as- 
suming the superposition of two Gaussian velocity distributions. 

This result was also co nfirmed over a w ider height interval, 
0.1 kpc < |z| < 3 kpc, by iLee et al. I (1201 ll) . who analysed the 
kinematics of ~ 17 500 G-type dwarfs from which they selected 
a subsample of bona fide thick disk stars by means of [Fe/H] 
and [cc/Fe] provided by the spectroscopic SDSS sample. The 
fact that othe r studie s based on the ph otometri c SPSS sample 
dlvezic et al. I (l2008l) . iBond et al. I (l2010l) . and lLoebman etakl 
(1201 ll) ) did not reveal any significant rotation-metallicity corre- 
lation seems to depend on the larger errors affecti ng the chemical 
abundances derived from the ugriz passbands (lArnadottir et al. I 



d2010h . ILeeet aTI d20TlK For these reasons, i t is of particu- 
lar importance to cite lKordonati s~et al. I (1201 ll) . who recently 
performed a spectroscopic survey of 700 thick disk stars us- 
ing fully independent VLT/FLAMES observations and found the 
same correlation within the interval 1 kpc < \z\ < 4 kpc from 
the plane, where the thick disk population dominates. Finally, 
a positive cor r elation was present in the sample discussed in 
iGratton et al. I d2003l) . who analysed 150 subdwarfs and early 
subgiants in the solar neighbourhood with homogeneous spec- 
troscopic abundances and accurate kinematics derived from the 
Hipparcos catalogue. 

The discussion above clearly shows that part of the burden is 
still on the observational side and therefore on the decisive con- 
tribution that Gaia' s parallaxes will bring to this t o pic, al ong the 
lines discussed in | Freem an. & Bland-Hawt horn I d2002l) and in 



iRe FiorentirTI (120121) for the MW halo. 



3. Chemical distribution 

We computed the expected gradient in [Fe/H] along the disk of 
the MW at various cosmic epochs. To do t hat we ran a detailed 
chem ical evolution model as described in Spitoni & Mat teucci I 
(1201 ll) . which is an upgraded version of the original two-infall 
model of Chiappini et al.(1997;2001), where details can be 
found. This model follows the evolution of 3 1 chemical species 
including H, He, D, C,N,0, ar-elements, Fe and Fe-peak ele- 
ments. Stellar nucleosynthesis is taken into account in detail 
and includes the contributions from low - and intermediate - 
mass stars (0.8-8MO), massive stars (M > 8M ) that die as 
supernovae (SNe) of Type II and Type Ib/c, and Type la SNe 
(white dwarfs in binary systems). One of the main assumptions 
of the model is that the thin disk of the Galaxy forms inside- 
out, namely the inner disk regions are assembled first by gas 
accretion of extragalactic origin (primordial chemical composi- 
tion) and then the outer disk regions form on timescales that vary 
from 2 Gyr in the inner disk (2 kpc) to 10 Gyr and more at larger 
galactocentric distances (<14 kpc). This ensures the formation 
of a gradient especially if coupled with a threshold in the gas 
density regulating star formation. Radial flows were not consid- 
ered [J 

The model adopted in this paper is Model S2IT of 
ISpitoni & Matteucci I (1201 ll) ( their Table 1): it is a two-infall 
model similar to tha t presented in Chiappini et al. (2001) and 
ICescutti et al. (2007). This particular model was chosen because 
it accurately reproduces the most recent estimate of t he present 
time abundance gradients bv lLuck & Lambert] (1201 ll) along the 
disk, who found A[Fe/H]/Afi = -0.062 + 0.002 dex kpc 1 . Our 
predicted present time gradient for Fe is A[Fe/H] / AR = -0.052 
dex kpc -1 in a galactocentric range 4-14 kpc. Given the uncer- 
tainties still existing in deriving chemical abundances, and in 
particular in the adoption of different model atmospheres, and 
considering the small difference in the galactocentric dis tance 
range adopted by us compared to ILuck & Lambert I (1201 ll) . we 
consider our model predictions to be in reasonable agreement 
wit h the observations. 

I Spitoni & Matteuccil (1201 ll) considered the gradient at 2 Gyr 
since the beginning of star formation and the one at 8 Gyr. While 
the gradient at 8 Gyr shows a clear increase of [Fe/H] with de- 
creasing galactocentric distance, the gradient at 2 Gyr shows 
an increase from the outer regions towards to 10 kpc, where it 



1 Including radial flows is necessary to generate abundance gradi- 
ents only if no inside-out and no threshold are assumed (Spitoni & 
Matteucci, 2011) 
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Table 1. Halos properties (Navarro, Frenk and White profile) 
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Fig. 1. Initial radial chemical distribution. The crosses show the model 
of Spitoni & Matteucci (2011), while the dots, which represent the 
chemical distribution of the N-body particles at t = Gyr, include an 
intrinsic dispersion cr lFc/H j =0.1 dex. 



reaches a peak, and then a decrease from R< 10 kpc, as shown in 
Figure 1 . This apparently strange behaviour was already present 
in the model of Chiappini et al. (2001) and can be explained 
on the basis of the inside-out disk formation, in particular, by 
the fact that at early epochs the efficiency of chemical enrich- 
ment in the inner regions of the disk is low due to the very 
large amount of infalling primordial gas. Then, at later epochs, 
while the star formation rate is still much higher than in the outer 
disk regions, the infall of primordial gas is much stronger in the 
outer disk. This produces the early flat gradient and the steepen- 
ing of the gradient in time. We adopted a metallicity dispersion, 
C[Fe/H] =0.1 dex, corresponding to the typical precision of the 
sp ectroscopic metallicity mea surements for SDSS, as estimated 
bv lAllende Prieto et al.Td2008h . 



4. Simulations 

We used a dark matter (DM) halo containing a rotating stellar 
exponenti al disk. The DM h alo has a Navarro, Frenk and White 
(NFW) (iNavarro et al.l[l9 97) radial density profile and structural 
parameters as indicated in Table 1 . These parameters are appro- 
priate for an MW-like DM halo at redshift z = 0. 

We used MW parameters at z=0 to initialise our model: since 
we want to compare gradients observed at present time with 
those generated by dynamical evolution, we need to have the 
correct physical parameters of today. Of course, this introduces 
a discrepancy, because several Gyr ago the Galaxy was differ- 
ent. However, we believe that our model, even if it is a simple 
nondissipative one, is adequate to separate effects of dynamics 
from chemical evolution. The DM particles of the nonrotating 
halo have velocities g iven by the local equilibrium approxima- 
tion dHern quist 1993). 

A stellar disk, whose mass is set to the value estimated for the 
MW, is embedded into our halo; the other features of the sim- 
ulated disk are described in Table 2. The position of each disk 
particl e is obtained by using the rejection method in lPress et al. I 
(1986); the disk is in gravitational equilibrium with the DM halo. 
We also simulated a comparison case consisting of the same disk 
but endowed with a bulge with a Hernquist radial density profile 
and mass of 1.8 x 10 M Q (see Table 3). The bulgeless disk de- 
cays in a barred configuration, while the disk with a bulge is sta- 
bilised against bar formation. We ran our two simulations using 



DM 


M vir 


Ryir 


r 




N 


e M DM 


Halo 


10 12 


258 


7.40 


336 


10' 


0.11 1.07 x 10 5 



Notes. M rir : Halo's virial mass in M ; R vir : virial radius in kpc. C,,„ : 
NFW concentration parameter. R max : maximum radius. N: total 
number of Halo particles, e: softening length in kpc. M DM : mass 
of DM particle in M G . 



Table 2. Properties of the disk 



Stars 


M, 


M slm 


h d 


Zd 


N 


Q 


Disk 


5.6-10 1 " 


7.47- 10 3 


3.5 


0.7 


7.5-10 6 


0.044 2 



Notes. M„: disk mass in M Q . M star : mass of the disk particle in M Q . 
hd'- disk scale length in kpc. Zd'- initial disk thickness. N: number of 
particles, e: softening length in kpc. Q: Toomre parameter. 

Table 3. Properties of the bulge (Hernquist profile) 



Stars 


M, 


N 


e a M h 


Bulge 


1.85-10 1 " 


2.5- 10 6 


0.044 1.12 7.39- 10 J 



Notes. M,:mass in M Q . N: number of particles, e: softening length in 
kpc. a: Hernquist scale radius in kpc. M b : mass of the particle in 

M . 



the public parallel treecode GADGET2 dSpringel. V. 1120051) on 
the cluster matrix at the CASPUR (Consorzio Interuniversitario 
per le Applicazioni del Supercalcolo) consortium, Rome. Both 
systems were left to evolve for 10 Gyr. The snapshots of the 
particle configurations were taken with a timestep of 0.25 x 10 8 
years. To avoid the effect of relaxation of our initial conditions, 
we let the dynamical system to evolve for 1 Gyr before assign- 
ing our initial metallicitie;fl and we analysed the redistribution 
of the chemical properties over the following 9 Gyr. We as- 
signed the radial chemical distribution function defined in Sect. 
[3] which is shown in Figure Q] to what we considered our initial 
configuration. Each particle in the configuration is tagged with a 
[Fe/H] label according to this initial radial function. 



5. Results 

5.1. Evolution ofV<p vs. [Fe/H] 

The evolution of the rotation-metallicity correlation 
dV<j,/d[F&/H] in the barred disk model, evaluated for parti- 
cles in the range -0.55 < [Fe/H] < -0.35 and > 50 km s _1 , 
is shown as a dashed thick line in Figures [3] |4] and [5] This 
choice for the chemical and kinematical ranges allows one to 
consider the core of the metallicity distribution and exclude the 
low-velocity tail. These figures represent the distribution vs. 
[Fe/H] at different times within the solar annulus 8 kpc < Rq < 
10 kpc and for 1.5 kpc < \z\ < 2.0 kpc. The colour-coded zones 
shown in the right panels of Figures [3] @] and [5] represent the 
distribution of the average (R ) original radii of the particles as 
a function of the location on the plane (V0,[Fe/H]). 
Table |4] lists the values of the correlation at various times and 
the time evolution of the correlation is shown in Figure [2] 



2 From this point on any reference to time evolution will be given 
with regard to the time of injection of the chemical distribution in the 
simulation. Therefore / = corresponds to 1 Gyr from the beginning of 
the dynamical simulation. 
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Fig. 2. Time evolution of the coefficient dV^/d[Fe/H\ 



By definition, at t — Gyr all particles belong to the annulus 8 
kpc < Rq < 10 kpc and no rotation-metallicity correlation can 
be present because the [Fe/H] values are randomly drawn and 
assigned to the particles at that time. Because of stellar motions, 
at t — 0.1 Gyr and in the 8-10 kpc annulus, we found many 
particles coming from the outer and inner regions, although no 
rotation-metallicity correlation is present as yet. 

A mild correlation appears between t = 0.2 and t = 0.3 Gyr, 
which quickly increases up to 82 + 6 km s _1 dex~' at t — 0.4 Gyr 
(Figures 2 and 3). We note that this correlation is produced by 
the "immigration" of low -rotating and metal-poor particles from 
the inner regions, which populate the bottom-left corner of the 
right panels (as indicated by the dotted white iso-contours) and 
tilt the iso-density contours in the left panels. As shown in Figure 
|U the rotation-metallicity distribution is stable over several Gyr, 
until about 6 Gyr. As reported in Table|4] the number of particles 
is also nearly constant at » 30000 in the 8-10 kpc annulus and 
1.5 kpc < \z\ < 2.0 kpc. 

After 6 Gyr, we observe a significant broadening of the 
distribution due to the arrival of many particles from the inner 
regions of the disk that populate the slowly rotating tail of the 
velocity distribution (V^ ^ 100 km s _1 ). We note, however, that 
a smaller and decreasing rotation-metallicity correlation persists 
up to t — 8 Gyr (eft. Table 3]). Notice also that the number of 
particles increases dramatically at 7 Gyr, reaching ~ 100000 
at t = 8 Gyr, with a corresponding dramatic change in the az- 
imuthal velocity distribution, as shown in the last two panels of 
Figure [6] 

This figure helps following the migration history of the parti- 
cles of interest by showing in the left panels the histograms of 
their original radial distributions at the same evolutionary times 
as those used for Figures [3] [4] and [5] Clearly the peak of the dis- 
tribution moves with the evolution of the disk outside of the solar 
annulus towards smaller initial radii. 

Finally, it is worth to remark on the intriguing feature of the 
initial sharp peak in Figure[2] We estimated an average timescale 
of « 0.2 Gy for the stars involved to oscillate around their guid- 
ing centres, a timescale compatible with the time needed for the 
formation of this initial spike in the evolution of the rotation- 
metallicity correlation. A deeper investigation on the role of the 
orbital features and on the dynamical orbital evolution will be 
the subject of a forthcoming study. 




Ro (kpc) 

Fig. 7. Upper panel: histogram of the original radial distributions of the 
subset of particles found within 8 kpc < R < 10 kpc after 5 Gyr. Lower 
panel: histogram of the original radial distributions of the same particles 
as in the upper panel, but with 1.5 kpc < |z| < 2.0 kpc. 



5.2. Barred vs unbarred disks, impact of metalliticy 
dispersion, and impact of relaxation 

Figure [7] upper panel, shows the original radial distribution of 
only those particles that after 5 Gyr are found within the solar 
annulus (8 kpc< R < 10 kpc), whereas the lower panel shows 
the subset of these particles that have 1.5 kpc < \z\ < 2.0 kpc. 
Here we can see how the outward radial migration is stronger in 
the barred disk: there, the bar i s dri ving the migration since, ac - 
cording to lFriedli etaD (1 19941) and iMinchev & Famevl d2010), 
the presence of nonaxisymmetric features in the disk causes ra- 
dial migration. On the other hand, it appears that a diffusion is 
also produced by the coarse-grained gravitational field present in 
our simulations, as indicated by the results of our nonbarred disk 
case. A fraction of this diffusion could be numerical; however, it 
might mimic th e diffusion generated by large molecular c louds 
in real galaxies dWielen I dl977h . |jenkins & Binnevl dl990l) ). 
In the barred-disk case the bulk of the particles at R = 8-10 kpc 
arrives from the inner region of the disk (Figure[7] upper panel). 
Because these particles originated in internal regions affected by 
a stronger vertical potential, they can extend to higher heights 
off the plane (Figure [7], lower panel). Th ese results are consis- 
tent with previous studies carried out by Roskaret al. (2008) 
and Loebman et al. 

In Figure [8] we show the histograms describing the distribution 
of the differences L/,„ - L,„, between the initial and final, after 
5 Gyr, angular momenta of the selected particles, in the solar 
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Fig. 3. Barred-disk model. Distribution of the azimuthal velocity and metallicity after an evolution of 0, 0. 1, 0.4 Gyr. Left panels: colour contour 
plot showing the rotation-metallicity distribution, vs. [Fe/H], of the particles in the solar annulus, 8 kpc < R < 10 kpc and 1.5 kpc < |z| < 2.0 
kpc. The dashed white line indicates the linear fit in the central region of the distribution, including all particles with -0.55 < [Fe/H] < -0.35 and 
V(j > 50 km s . Right panels: the colour contour plot represents the distribution of the mean original radius, (Rq), as a function of ([Fe/H],V<j) of 
the same particles shown in the left plot. Here, the dotted white lines mark the iso-contours at {R )=6 kpc, 8 kpc, 10 kpc, and 12 kpc. The black 
contours indicate the iso-density levels of the ([Fe/H],^) distribution shown in the left panel. 
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Fig. 4. Distribution of the azimuthal velocity and metallicity after an evolution of 1, 3, 5 Gyr. 



annulus 8-10 Kpc, for the barred and unbarred disks. In the gular momentum exchanges are more important for the barred 

top histograms, we present the angular momenta differences for case. The less evident skewness of the histogram for the barred 

all particles in the annulus, whereas the lower panels show the case and particles at higher \z\ reflects the more symmetric dis- 

same differences for these particles with heights above the plane tribution shown by the red curve in the lower panel of Figure 

in the interval 1.5 kpc < |z| <2 kpc. These plots show how an- [7] 
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Fig. 5. Distribution of the azimuthal velocity and metallicity after an evolution of 6, 7, 8 Gyr. 



To test the effect of metallicity dispersion on our results, we 
repeated our numerical experiment by plugging the same metal- 
licity profile, but with different initial metallicity dispersions. As 
expected, the rotation-metallicity correlation is strongly depen- 
dent on the assumed initial metallicity dispersion. Indeed, at t=5 



Gyr its value becomes 21+8 km s dex -1 if we run our model 
with an initial metallicity dispersion of 0.2 dex, while it becomes 
157 + 6 km s dex -1 with an initial metallicity dispersion set to 
0.05 dex. 

We note that the impact on the rotation-metallicity correlation of 
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Fig. 6. Left panels: Normalised histograms of the original radii of the subset of particles found within 8 kpc < R < 10 kpc that have 
1.5 kpc < \z\ < 2 kpc at different evolutionary times. Right panels: Normalised distributions of the azimuthal velocity at the same 
times. 



a higher initial metallicity dispersion (a reduction in the corre- 
lation) appears to be similar to that caused by churning, which 
might be dominating after 6 Gyr (see Figure[2]l. This will be in- 
vestigated in our forthcoming study. 

Finally, to address the question of the possible role of par- 
tial/incomplete dynamical relaxation, we re-assigned the initial 
metallicities after 2 Gyr of evolution and compared the results 
to the 1-Gyr case. No significant difference was found, thus con- 
firming that our results are not the spurious effect of relaxation 
of the adopted initial conditions. 

5.3. Comparison to the data 

The results of Figure 2 testify to a positive rotation-metallicity 
correlation that persists at ss 60 km s _1 dex -1 for several Gyr, 
to be compared wi th t he 40-50 km s" 1 dex -1 measured by 
ISpagna et al.l d20ld) and lLee et al. I d20TTh . 
Although we did not attempt any rigorous fit to the data, which 
would not be appropriate given the limitations of our galaxy 
model, we emphasise the consistency of the kinematics pro- 
duced by our dynamical simulations with what was unveiled by 
the data. 

Figures l9l and [TOl compare the velocity distribution, vs. Vr, 
of the A-body particles in the peripheral z-layers to the stellar 
dis tribution derived from the SDSS - GSC-II catalogue studied 
by|Spagna et al. (2010). The theoretical distribution shows par- 



ticles from the adopted barred-disk model evolved for 5 Gyr and 
selected within the solar annulus at 1.5 kpc < \z\ < 2.0 kpc. We 
note that, in spite of the model limitations, the two distributions 
appear to be quite similar, the moments of the simulated distribu- 
tion ((V^) - 142 km s , o>0 - 59 km s , o> (( ^ 70 km s _1 ) are 
quite close to the values of the observed distribution ((V^) - 159 
km s , o>0 a 56 km s _1 , cry R - 62 km s . 
In addition, we note an excess of slowly rotating particles (V^ ^ 
100 km s _1 ) with respect to the observations; these are parti- 
cles coming from the innermost regions of the disk, and this ex- 
cess seems to depend on the bar generated in the simulated disk, 
which is stronger than the MW bar and likely more efficient in 
displacing particles. 

6. Discussion and conclusions 

Nondissipative mechanisms inside disk s due to an irreg ular grav- 
itational field as already theorised by IWielen I (Il977h are very 
efficient in promoting stellar orbital diffusion. Nonaxisymmetric 
structures s uch as spiral arms and bars prom ote significant radial 
migra ti on (dSellwood & Binnev 1120 02) and ( Minchev & Famevl 
2010)). ISchonrich & Binnevl (12009a) discussed in detail the 
churning and blurring processes that are responsible for the ra- 
dial mixing. Through the comparison between barred and un- 
barred disk models, our simulations confirm the orbital diffusion 
in both cases but a more efficient radial mixing in the presence 
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Fig. 8. Normalised distributions of the differences Ly,„ - L,„, between initial and final total angular momenta (after 5 Gyr) of the particles inside 
the solar annulus for the barred and unbarred disk (the curves represent an interpolation of the histogram). 



of a bar. These dynamical processes, which naturally produce 
secular migration of the stars from the disk central regions to- 
wards larger radii and higher vertical distances from the plane, 
represent a plausible mechanism to explain the origin of the MW 
thick disk. 

We recall that our simulations do not include any gaseous com- 
ponent or star formation. We simply investigated the evolution 
of a system formed by collisionless particles representing the 5- 
10 years old parent population of the thick disk component we 
observe today. Through this approach, we dramatically reduced 
the number of degrees of freedom of the model, so that we were 
able to investigate the details of the dynamical processes and 
their effects on the thick disk chemical properties in situ, i. e., at 
higher heights from the plane, where we assume that no signifi- 
cant contamination due the younger thin disk stars is present. In 
the lig ht of our results, the thick disk rotatio n-metallicity corre- 
lation (Spagna et al. 1 120 lot iLee et al. 11201 lb seems to represent 
an important signature of the disk evolution. This feature can in- 
deed be explained by the role of the dynamical migration in our 
iV-body simulations, once we assume an initial radial chemical 
gradient as that suggested b y the chemical evolution models of 
Spitoni & Matteuccy d201 lh . which prescribes a positive slope 
for the inner early disk, R ^ 10 kpc, combined with the usual 
decreasing slope in the outer disk. The crucial role of a posi- 
tive inner slope for the primordial chemical distribution to pro- 
duce a positive rotation-metallicity correlation was presented in 



ICurir et aT. 1 (120121) . where a distribution consisting of two sim- 
ple linear functions with positive slope for inner radii up to 6 kpc 
and a negative one outwards was plugged in the early (N-body) 
disk.0 

The sce nario above i s consi stent with the recent results pub- 
lished by ICresci et"aL l d20Toh . who studied a sample of distant 
galaxies at redshift z > 3 and found evidence for an inverse 
metallicity gradient, possibly produced by the accretion of pri- 
mordial gas, which diluted the abundance of elements heavier 
than helium in the centre of the galaxies. Thu s, as the negative 
correlation, -22 km s dex , of the thin disk dLee et al. 1 l201lh 
results from the current - monotonically decreasing - radial gra- 
dient, the positive rotation-metallicity correlation of about 40-50 
km s _1 dex 1 , shown by the ancient thick disk population, can be 
explained by the relic signature of the inverse chemical gradient 
of the early ISM from which these stars formed 8-12 Gyr ago. 

Acknowledgements. A.C. and P.R.F. acknowledge the financial support of the 
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Mission - The Italian Participation to DPAC). A.C. and G.M. acknowledge the 
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an Italian Network for Computational Cosmology". The simulations were car- 
ried out at CASPUR, with CPU time assigned through the Standard HPC grant 
2010 : "Secular processes in disk galaxies and the formation of the thick disk". 



3 This proves that the radial location of the peak of the metallicity 
distribution adopted in this paper has no relevance for the formation of 
the rotation-metallicity correlation. 
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Table 4. Evolution of the rotation-metallicity correlation 



300 



200 



I 100 r 




100 



-300 -200 -100 100 200 300 

V r (km s" 1 ) 

Fig. 9. Simulated distribution of the azimuthal velocity vs. radial ve- 
locity based on particles (8 kpc< R < 10 kpc and 1.5 kpc < |z| < 2.0 
kpc) from the barred disk model at f = 5 Gyr. Here, blue and red sym- 
bols mark the particles with simulated [Fe/H]> -0.5 and [Fe/H]< -0.5, 
respectively. 
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Fig. 10. Distribution of the azimuthal velocity vs. radial velocity based 
on 3507 stars of the SDSS - GSC-II catalog (1.5 kpc < |z| < 2.0 kpc 
and -1.0 < [Fe/H] < -0.3). Here, blue and red symbols mark the stars 
with measured [Fe/H]> -0.6 and [Fe/H]< -0.6, respectively. 
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